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determine a distribution over avalanche ensemble with negative temperature. 
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noise of the energy plays a crucial role. This distribution is shown to be 
solution of both fractional Fokker-Planck equation and nonlinear one. As a 
result, we obtain new relations between exponent of the size distribution, frac- 
tal dimension of phase space, characteristic exponent of multiplicative noise, 
number of governing equations, dynamical exponents and nonextensivity pa- 
rameter. 
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I. INTRODUCTION 



In recent years considerable study has been given to the theory of self-organized criti- 
cality (SOC) that explains spontaneous (avalanche-type) dynamics, unlike the typical phase 
transitions that occur only when a control parameter is driven to a critical value fll], Q. A 
main feature of the systems displaying SOC is their self- similarity that derives to a power-law 
distribution over avalanche sizes. Respectively, SOC models are mostly studied by making 
use of the scaling-type arguments supplemented with extensive computer simulations (see 
0]). On the contrary, in this paper we put forward an analytical approach, which is enable 
to describe in phenomenological manner both process of a single avalanche formation and 
behavior of whole avalanche ensemble. 

The SOC behavior appears in vast variety of systems such as a real sand pile (ensemble of 
grains of sand moving along increasingly tilted surface) 0] - 0]) intermittency in biological 
evolution |^, earthquakes and forest-fires, depinning transitions in random medium and 
so on (see 0). Among the above models the sandpile is the simplest and most widely 
studied as analytically ||lOl, so numerically [jl2| - In analytical direction, a variety 



of field theory approaches is worthwhile to notice. Among them, the field scheme |T5 
based on nonlinear diffusion equation has failed to account for the main feature of self- 
organizing systems - the self-consistent character of avalanche dynamics. The visible reason 
is that using one-parameter approach does not take into account a feedback between open 
subsystem and environment that are related to order and control parameters, respectively 
(see criticism in Refs. [Q, also). Much more substantial picture is given within two- 
parameter approaches - ? EH that use both fundamental fields: gauge ones related to 
hydrodynamical modes type of sand pile height and material fields as a number of moving 
sand grains (avalanche size). Making use of the mean- field approximation shows that the 
self-similar regime of the sand pile dynamics is relevant to subcritical behavior, where a 
characteristic time of the order parameter variation is much more than the same for the 
control parameter, and the latter follows the former adiabatically. Such type adiabatic 



behavior is inherent in usual regime of system evolving in course of phase transitions |T7 
and jammed motion of vehicles |T3, so that adiabatic approach will be put on basis of our 
consideration. 

Perfect treatment of the SOC has been achieved within three-parameter approach based 
on the Reggeon field theory that uses the density of active sites pa as the order parameter 



and conserved field of the energy density C, as the control parameter |]I9|, [^. Along this 
line, SOC regime appears as a result of competition between a rate of the energy input h > 
and a dissipation rate e. It is appeared that the system under consideration behaves in quite 
different manner in case of fixed energy, when h = e = and total energy is conserved, and 
for driven sandpile, when h — >■ O"*", e — >■ at stationarity condition e > h. The first case shows 



TP| , to be reduced to the picture of supercritical regime, where nonhomogeneity of initial 



distribution of energy arrives at a non-Markovian term and space-dependent parameters of 
theory. At dimensions above critical value dc = 4, this case is found to be identical to the 
simplest Landau picture with ~ C ~ Cc in active stationary state {( > Cc) and Pa = 
in absorbing configuration {( < (c)- Fundamentally different picture appears in the case 
of driven sandpile, where due to external input /i — > O"*" the energy density is lost as an 
independent field being reduced to the critical value Cc- Here, the average magnitude of 
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the density of active sites is equal {pa) = h/e, so that the susceptibihty x = (dpa/dh) is 
turned out to be x = Respectively, a response function behaves at large distances r as 
x(r) ~ r'^-'^Q-^/^^ where c? is a space dimension, ^ ~ e~^/^ is a correlation length being scale of 
a system size L ~ e~^/^. It is remarkable that such mean- field-type behavior is caused solely 



by stationary condition and translational invariance [Q. A set of basic critical exponents 
governing scaling avalanche formation reads [|l^]: /5 = 7 = 5 = 1, /i = 2, z/ = l/2 and 77 = 0. 
On the other hand, scaling relations accompanied by equality of the susceptibility and mean 
size of avalanche yield the following expressions 

for the exponents of the avalanche size distribution 

P^s,e) = s-^V{x)- x = s/s,, s,~e-l/^ (2) 

where a critical size Sc is connected with system size L ~ ^ and a characteristic time 
tc ~ as Sc ~ ~ t^^^ (exponents D = p/a and z are fractal dimension and dynamical 



exponent related to a critical avalanche, respectively). According to Ref. [0, the mean- field 
magnitudes of the above exponents are as follows: r = 3/2, a = 1/2, D = 4 and z = 2. 

Along the standard approach |^ , we will use as gauge, so material fields: being reduced 
to velocity components and sand pile slope, the formers are considered at study of a single 
avalanche formation, whereas the latters are reduced to the number of moving sand grains 
at examination of a distribution over avalanche sizes. Section |I| contains the self-consistent 
theory of the flux steady-state developed along the first direction. It enables us to treat 
the problem of a single avalanche formation on the basis of the unified analytical approach 
being relevant to the case of fixed energy in Refs. ^U\. Section |T| deals with accounting 



additive noises of the sand velocity components and sand pile slope. The noise intensities 
increase is shown to make possible avalanche emergence even in non-driven systems. By 
this, the control parameter noise plays a crucial role. Such type fluctuational regime cor- 
responds to the case h ^ 0^ |^0|, where a distribution of order parameter is appeared in 
an algebraic form with integer exponent. In order to be not restricted by such particular 
case, we introduce unified Lorenz system with fractional feedback in Section Using this 
supposition allows us to describe subcritical regime of the avalanche formation in natural 
manner. The above generalization is put forward basis of Section devoted to consider- 



ation of avalanche ensemble. Following famous Edwards paradigm [E2], PB[, an effective 



scheme addressed to nonextensive thermodynamics ||2^ is proposed to be determined by a 
time-dependent distribution over energies of moving sand grains. To generalize the Edwards 
scheme to nonstationary nonextensive system, we use the fractional Lorenz system, where 
the avalanche size plays a role of the order parameter, nonextensive complexity is reduced to 
the conjugate field and nonconserved energy of the moving grains is the control parameter. 
Within the framework of such approach, the phase diagram is built to define the different 
domains of system behavior as a function of noise intensities of above values. As a result, we 
arrive at natural conclusion that the power-law distribution inherent in the SOC regime 
is caused by noise of the energy. 

Section |VI| shows that this distribution is the solution of both nonlinear Fokker-Planck 



equation, which appears at studying nonextensive systems p4[, and fractional Fokker-Planck 
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equation inherent in Levy-type processes characterized by dynamical exponent z |2^. As 
a result, we obtain new relations between the exponent r of the distribution (j^), fractal 
dimension D of phase space, characteristic exponent of multiplicative noise, number of gov- 
erning equations being needed to present self-consistent behavior in SOC regime, dynamical 
exponent z and Tsallis nonextensivity parameter q. 

Section |VIII| is Appendix containing basic properties of fractional integral and derivative. 



as well as Jackson derivative. 



II. NOISELESS AVALANCHE FORMATION 

Within framework of the simplest model of a real sand pile, a dependence y = y{t, x) 
defines its surface at given instant of time t. Locally the flow of sand can be described 
in terms of three quantities: the horizontal and vertical components of the sand velocity, 
X = dx/dt, y = dy/dt, and the surface slope y' = dy/dx. The key point of our approach 
is that the above degrees of freedom are assumed to be of dissipative type, so that, when 
they are not coupled, their relaxation to the steady state is governed by the Debye-type 
equations: 

dx X , , 



# ^ y_ 

dt rr 

dy' y'o - y' 



dt Ts 



(4) 



(5) 



where r^, t^"* and Ts are the relaxation times of the velocity components and the slope, 
respectively. Eqs. @ - (j^) imply that the sand is at rest in the stationary state, x = ij = Q 
and the equilibrium slope y' = y'o plays the role of a control parameter. 

Since the motion of sand grain along different directions is not independent, Eq. 
should be changed by adding the term / = y/7 due to liquid friction force along the y-axis 
(7 being the kinetic coefficient). Then, we have 

TxX = -X + a~^y, (6) 

where a = Note that, owing to the diffusion equation y = Dy" {D is the diffusion 

coefficient), the friction force appears to be driven by the curvature of the sand pile surface: 

/ = {D/i)y". (7) 

On the other hand, when x = (stationary state), solution of Eq. @ defines the tangent 
line y = ax + const, so that the friction force / = t~^x is proportional to the horizontal 
component of the sand velocity. Taking into consideration the relation and using the 
chain rule dy'/dt = y' + y"x, from Eq. one obtains the equation of motion for the slope: 
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-^sy = {y'o - y) - i^s/D) yx. 



(8) 



Following the same line, the equation for the vertical component of the velocity can be 
deduced 

r„S = -, + ^,'i. U±^(l + fJpi. (9) 

Note the higher order terms are disregarded in Eq. (|^) and the renormalized relaxation time 
Ty depending on the stationary slope y'^ is introduced. 

Equations (|), (|) obtained constitute the basis for self-consistent description of 
the sand flow on the surface with the slope y' driven by the control parameter y'^. The 
distinguishing feature of these equations is that nonlinear terms, that enter Eqs. (|^), (^, 
are of opposite signs, while Eq. is linear. Physically, the latter means that on the early 
stage the avalanche begins motion along the tangent y = ax + const. The negative sign of 
the last term in Eq. (|]) can be regarded as a manifestation of Le Chatelier principle, i.e., 
since the slope increase results in the formation of an avalanche, the velocity components x 
and ij tend to impede the growth of the slope. The positive feedback of x and y' on ij in 
Eq. @ plays a fundamental part in the problem. As we shall see later, it is precisely the 
reason behind the self-organization that brings about the avalanche generation. 

After suitable rescaling, Eqs. @, (H), can be rewritten in the form of the well-known 
Lorenz system: 

ii = —u + w, (10) 
ei) = —V + uS, (11) 
SS = {So -S)- uv, (12) 

where u={Ty/Txy^'^{Ts/Dy^'^x, v={Ty / r^^y/'^ {ts / DY^'^ij / a and S={Ty/Tx)y'/a are the dimen- 
sionless velocity components and the slope, respectively; e = Ty/T^, 5 = ts/t^ and the dot 
now stands for the derivative with respect to the dimensionless time I/tx- In the general 
case, the system (|I0|)-(|1^) cannot be solved analytically, but in the simplest case, where 
e -C land 5 -C 1, the vertical velocity v and the slope S can be eliminated by making use 
the adiabatic approximation that implies neglecting of the left hand side of Eqs. (pT]), (|12|). 
As a result, the dependencies of S and v on the horizontal velocity u are given by 



c "^0 SqU 



Note that, under u is within the physically meaningful range between and 1, the slope 5* is 
a monotonically decreasing function of u, whereas the velocity v increases with u {aX u > 1 
we have dv/du < and the flow of the sand becomes unstable). 

Substitution of the second Eq. (|13D into Eq. (p!0D yields the Landau-Khalatnikov equation 



du 

with the kinetic energy 



« = -?^ (14) 
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E = -u' 



u 



(15) 



For 5*0 < 1, the w-dependence of E is monotonically increasing and the only stationary 
value of u equals zero, uq = 0, so that there are no avalanches in this case. Obviously, such 



steady-state relevant to absorbing configuration studied in Ref. [Q. If the slope 5*0 exceeds 
the critical value, Sc = 1, the kinetic energy assumes the minimum with nonzero steady 
state velocity components Ue = = {Sq — 1)^/^ and the slope 3^ = 1- 

The above scenario represents supercritical regime of the avalanche formation and ad- 
dresses to the second-order phase transition |T^. The latter can be easily seen from the 



expansion of the kinetic energy (|T^) in power series of <^ 1. So the critical exponents 



7, S, V are identical to those obtained within the framework of the mean-field theory [p!9 
However, the magnitude /? = 1/2 is twice as little because our order parameter (being the 
velocity) is not reduced to the same (the number of active sites) in theory [^. 

The drawback of the outlined approach is that it fails to account for the subcritical regime 
of the self-organization that is the reason for the appearance of avalanches and analogous to 
the first-order phase transition, rather than the second-order one. So one has to modify the 
above theory by assuming that the effective relaxation time t,j.[u) increases with velocity u 
from value '^^.(l + m)~^, m > to r^. 11^. The simplest two-parameter approximation is 



1 + 



1 + (u/uo)^ 



(16) 



where < mq < 1. The expression for the kinetic energy ([15D then changes by adding the 
term 



AS 



Un In 



1 + 



u 



(17) 



and the stationary values of u are following: 



Mo 



ZUqq 



Moo |1 =F [1 +^0^00 i^o 

(5*0 - 1) + Scul, Sc 



Sr 



1/2 



m. 



1/2 



(18) 



The upper sign in the right hand side of Eq. (|T^) meets the value of the unstable state m™, 
where the kinetic energy E + AE has the maximum, the lower one corresponds to the stable 
state Ue- The corresponding value of the stationary slope 



1 + ^2 



00 



00) 



0) ^0 



smoothly increases from the value 



Sr, 



at the parameter = Sco with 



1 +Uo\Jm/{l -ul) 



S 



cO 



(l - ul) 



(19) 



(20) 



(21) 
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to the marginal value 5'c = 1 + m at 5*0 = Sc- The 5'o-dependencies of Me, Um, and 5*6 are 
presented in Fig. 1. As is shown in Fig. la, under the adiabatic condition ts <^ is met 
and the parameter 5*0 slowly increases being below Sc {Sq < Sc), no avalanches can form. At 
the point 5*0 = Sc the velocity Uc jumps upward to the value V^Uqq and its further smooth 
increase is determined by Eq. (p^. If the parameter 5*0 then goes downward the velocity Ue 
continuously decreases to the point, where 5*0 = ^coandug = uqq. At this point the velocity 
jump-like goes down to zero. Referring to Fig. lb, the stationary slope Se shows a linear 
increase from to Sc with the parameter 5*0 being in the same interval and, after the jump 
down to the value (1 — Mq)"^ '^'o = Sc, Se smoothly decays to 1 at 5*0 3> Sc- Under the 
parameter So then decreases from Sc down to Sco, the slope grows. When the point (pT]) is 
reached, the avalanche stops, so that the slope undergoes the jump from the value ( pO|) up 
to the one defined by Eq. pl|). For 5*0 < Sco again the parameter Se does not differ from 
Sq. Note that this subcritical regime is realized provided the parameter m, that enters the 
dispersion law (|16|), is greater than 

m^in = (22) 

Clearly, according to the picture described, the avalanche generation is characterized by 
the well pronounced hysteresis, when the grains of sand initially being at rest begin to move 
downhill only if the slope of the surface exceeds its limiting value Sc = 1 + m, whereas the 
slope Sco needed to stop the avalanche is less than Sc (see Eqs. (pUD, (PH)). This is the case 
in the limit ts/t^ and the hysteresis loop shrinks with the growth of the adiabaticity 
parameter 6 = ts/t^. In addition to the smallness of S, the adiabatic approximation implies 
the ratio Ty/Tx = e is also small. In contrast to the former, the latter does not seem to be 
realistic for the system under consideration, where, in general, Ty ^ t^. Thus, it is of interest 
to study to what extent the finite value of e could change the results. 

Owing to the condition 5 -C 1, Eq. (0) is still algebraic and S can be expressed in terms 
of u and v. As a result, we derive the system of two nonlinear differential equations that 



can be studied by the phase portrait method |T^. The phase portraits for various values of 



e are displayed in Fig. 2, where the node point O represents the stationary state and the 
saddle point S is related to the maximum of the kinetic energy. As is seen from the figure, 
independently of e, there is the universal section that attracts all phase trajectories and its 
structure is appeared to be almost insensitive to changes in e. Analysis of time dependencies 
v{t) and u{t) reveals that the velocity components slow down appreciably on this section in 
comparison to the rest parts of trajectories that are almost rectilinear (it is not difficult to 
see that this effect is caused by the smallness of the parameter 6) . Since the most of time the 
system is in vicinity of this universal section, we arrive at the conclusion that finite values 
of e do not affect qualitatively the above results obtained in the adiabatic approximation. 



III. NOISE INFLUENCE ON AVALANCHE FORMATION 

We focus now on studying the affect of additive noises of the velocity components u, 
V, and the slope 5*. With this aim, we should add to right hand parts of Eqs. ([I0| ) - (|12D 
the stochastic terms I^'^C,, l]^'^^, ^s^"^^, respectively [here the noise intensities /„, 1^, Is are 
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measured in units {Tx/Ty){D /ts), a'^ {t^ / Ty) {D / ts) , a^ijx/Ty)^ correspondingly, and ^{t) is 
5-correlated stochastic function] Then, within the adiabatic approximation, equations 
(|n|), (0) could be reduced to the time-dependencies 

v{t)=v + vm, S{t) = S + Sm; (23) 



v= Soud{u), V = \J Iv + Isu^ d{u), (24) 



2\-l 



S= Sod{u), S = y Is + IvU^ d{u), d{u) = {1 + u 

Here deterministic components are reduced to Eqs. whereas fiuctuational ones follow 
from the known property of additivity of variances of independent Gaussian random quanti- 
ties 1 5^. Thus, using slaving principle inherent in synergetics transforms noises of both 



vertical velocity component v and slope S, having be adiabatic initially, to multiplicative 
form. As a result, combination of Eqs. (0), (^1), (|2^) arrives at Langevin equation 

dE 



u = f{u) + ^I{u)m, / = (25) 

where force / is related to the energy E determined by Eq. (|T3|) and expression for effective 
noise intensity 

I{u)^I^+[h + Isu')d\u) (26) 

is obtained in accordance with above mentioned property of noise variances additivity. In 
order to avoid mistakes, one should notice that direct insertion of Eqs. (p3D, into (p!0[) 
results in appearance of stochastic addition 

i'j' + {iy'+i'J'u)d{u)]m, (27) 



whose squared amplitude is quite different from effective noise intensity (|26|). Moreover, con- 
trary to expressions (p4D, direct using adiabatic approximation in Eqs. (0), ([l^) reduces fiuc- 
tuational additions in Eqs. ( pSj) to the forms: v = (jy^ + ll^'^u)d{u), S = (/^^ — I^^'^u)d{u). 
The latter is obviously erroneous since the effective noise of the slope S disappears entirely 



for the horizontal velocity u = yls/Iv The reason of such contradiction is that Langevin 
equation, being stochastic in nature, does not permit usage of usual analysis methods (see 

To continue along standard direction, let us write the Fokker-Planck equation related to 



Langevin equation (|25D : 

'-^^l{-nu)Piu.t)^llHu)PM]]. (28) 

At steady state, which will be considered only, the probability distribution P{u, t) becomes 
time-independent function P{u) and usual condition, that expression in braces of right hand 
part in Eq. ( P^ is zero, arrives at stationary distribution 

P{u) = Z-^ exp{-U{u)}, (29) 
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where Z is a normalization constant and effective energy 



f{u') dE 
— -du, / = -— , 
[u') ou 



(30) 



is determined by the bare energy E, Eq. ( ]T5| ) and the noise intensity /(m), Eq. 
Combining these expressions, we might find the exphcit form of U{u), which is too cum- 
bersome to be reproduced here. The equation which defines the locations of the maxima of 
distribution function P(u) 



- Sqx^ - 2Isx + 4(/s - 7^,) = 0, x = l + u'^ 



(31) 



is much simpler. According to Eq. (0), they are insensitive to changes in the intensity of 
noise /„ of the velocity component u, but are determined by the value 5*0 of the sand pile 
slope and the intensities ly, Is of the noises of vertical velocity component v and slope S, 
which acquire the multiplicative character in Eq. (^). Hence, it can be put for simplicity 
= and Eqs. ([15|) , (|30D , (|26|) give the following expression for the effective energy: 



U{u) 



u 



+ (2 - 5o - i)u^+ 



(32) 



(1 - 2) (1 - 5o - i) Hi + u^) + Is ln[g(u) + igl{u)i i = I Jig. 



According to Eq. (|3TD, the effective energy (|32D has a minimum at m = if the driven 
slope Sq does not exceed the critical level 



Sc = l + 2Is 



4/„ 



(33) 



whose value increases at increasing intensity of noise of the sand pile slope, but decreases 
with one of the vertical velocity. Here, sand grains are at rest. In the simple case ly = 0, 
the avalanche creation is related to solutions 



5o - 3 + V(3 - Soy + 4(2So - 3 + 2/s 



(34) 



which are obtained from Eq. ( 
this solution has its minimum 



after elimination of the root = 0. The magnitude of 



{So - 3) - v/(5o + 7)(5o-l) 



(35) 



on the line defined by expression (l33| ) with ly = 0. At 5*0 < 4/3 the roots ±Uc are complex, 
at 5*0 = 4/3 they become zero, at 5*0 > 4/3 they are real and m+ = — m_. In this way, the 
tricritical point 



^0 = 4/3, 



1/6 



(36) 



addresses to the appearance of roots M-t 7^ of Eq. (|3lD (avalanche creation). If condition 
(PB|) is satisfied, the root u = corresponds to the minimum of the effective energy (B^) at 
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5*0 < 4/3, whereas at 5*0 > 4/3 this root corresponds to the maximum, and the roots u± - 
to symmetrical minima. 

Now, we find another condition of stabihty of roots u±. Setting the discriminant of 
Eq. ( ^11) equal to zero, we get the equations 



ls = 0, 

the second of which gives 



Is-Is 



27 



So 

y 



C2 
•^0 



+ 



Sn 



0, 



(37) 



21, 



± 



1/2 



2S'o 



(38) 



This equation defines a bell-shaped curve So{Is), which intersects the horizontal axis at the 
points Is = and Is = 27/2, and has a maximum at 



5*0 



(39) 



It is easy to see that for ly = this line touches the curve (|33D at point (pGl). 

Let us consider now the more general case of two multiplicative noises Iv,Is 7^ 0. In- 
troducing the parameter a = 1 — i, i = I^/Is and the renormalized variables / = Is/o?, 
So = So/ a, V? = {1 + v?) / a — \, at z < 1 we reproduce all above expressions with the general- 
ized energy U / 1 in Eq. (0). Thus, action of the noise of the vertical velocity component v is 
reduced to the renormalization of the extremum value of the horizontal one by the quantity 
(a~^ — 1)^/^. As a result, the region of divergence £t ^ becomes inaccessible. 

The condition of extremum of the generalized energy ( ^2t) splits into two equations, one 
of which is simply m = 0, and the other is given by Eq. (|3l|). As mentioned above, analysis 
of the latter indicates that the line of existence of the zero solution is defined by expression 
(0). The tricritical point has the coordinates 



5o 



:i - /„ 



6 



:i + 84 



(40) 



The phase diagram for the fixed intensities I^ is shown in Fig. 3. Here the curves 1, 2 define 
the thresholds of absolute loss of stability for the fiuxless and fiux steady-states, respectively. 
Above line 1 the system occurs in a stable flux state, below curve 2 it is in fluxless one, and 
between these lines the two-phase domain is realized. For < 1/4 situation is generally 
the same as in the simple case = (see Fig. 3a). At J^, > 1/4 the avalanches formation 
is possible even for small intensities Is of the slope noise (Fig. 3b). According to (|iO|), the 
tricritical point occurs on the Js'-axis at It, = 1, and for the noise intensity 1^ larger than 
the critical value = 2 the stable fluxless state disappears (see Fig. 3c). 

The above consideration shows that the dissipative dynamic of grains flow in a real sand 
pile can be represented within the framework of the Lorenz model, where the horizontal 
and vertical velocity components play roles of the order parameter and its conjugate fleld, 
respectively, and the sand pile slope is the control parameter. In Section the noiseless 
case is examined to show that an avalanche is created if the externally driven sand pile slope 
is larger than the critical magnitude 
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(41) 



In this sense, the systems with small values of the kinetic coefficient 7 and large relaxation 
times Tx, Ty of the velocity components are preferred. However, the sand flow appears here 
as a phase transition because the avalanche creation in the noiseless case is possible only 
due to the externally driven growth of the sand pile slope. 

Accounting the additive noises of the above degrees of freedom shows that the stochas- 
ticity influence is non-essential for the horizontal velocity component and is crucial for both 
the vertical one and the sand pile slope. The avalanche appears spontaneously if the dimen- 
sionless noise intensities are connected by linear relation 

Is = -l + 2h, (42) 



following from Eq. (^) at conditions a; = 1 (m = 0), 5*0 = 0. According to Eq. (^2|) 



in absence of the sand pile slope noise the avalanche is created if the intensity of vertical 
velocity component exceeds the value 

1 

/.o = J (43) 

corresponding to the point O in Fig. 4. Increase of both the vertical velocity and the sand 
pile slope noises causes avalanche formation if its intensities are bounded by the condition 
(1^). With further increase of these intensities above magnitudes 

hi = ^^, hi = -^, (44) 

the domain of the mixed state appears at the point T in Fig. 4. If the noise intensity of the 
vertical velocity is more than the larger value 

1,2 = 2 ^, (45) 

TxTyTs 

corresponding to the sand pile slope noise 1 32 = Gl'^/TxTy (the point C in Fig. 4), the fiuxless 
steady- state disappears at all. 

Physically, we have to take into consideration that the SOC regime is not relevant to 
fiux-type avalanche state itself, but rather to an intermittent regime of avalanche formation 
corresponding to the domains on the phase diagrams in Figs. 3, 4, where a mixture of both 
phases A and N (avalanche and non-avalanche) exists. According to above analysis, such 
an intermittent behavior may be realized within the region located upper the line (|42|) and 
outside the curve determined by the equation 



V 



V27 



(46) 



for the dimensionless values /„, Is- Corresponding phase diagram is depicted in Fig. 4 to 
show very non-trivial form (especially, within the domain I^i < I, < 1,2)- 
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IV. GENERALIZING SELF- SIMILARITY 



To proceed the consideration of the system behavior, let us examine exphcit form of the 
probabihty ( p9[ ) determined for different regimes by the effective energy (|30|). In the case 
Iu,Is ^ we obtain Gibbs-type distribution 

Piu) ^ + u'Y exp J /(«)(! + u'fdu^ , f{u) ^-u + S,u/{1 + u"), (47) 

being opposite to the power dependence inherent in self-similar systems. Contrary, at inter- 
mittent behavior, when /i, ^ Isi supercritical values of the slope noise intensity Is cause 
the following distribution form: 

Thus, the case <^ Is addresses to the power-law distribution being relevant to self- 

similar behavior. However, as corresponding consideration shows, obtained exponent is 
not reduced to integer 2, generally. 

To get rid off such a restriction, the multiplier u in nonlinear terms of Eqs. (|TU|) - (|T2p 
is supposed to be replaced by m", where an exponent < a < 1. With accounting the 
stochastic additions one obtains then the basic equations in dimensionless form 



U = -U + V + \J lu ^(t), 

ev = -v + u''S+^,i{t), (49) 
5S = {So- S)-u''v + JTs i{t). 



Thus, it is appeared that the agreement of the Lorenz self-organization scheme with SOC 
conception related to self-similar systems is achieved if one assumes that both positive and 
negative feedbacks are fractional in the nature. Within such a supposition, the adiabatic 
approximation e,5 <^1 arrives at the Langevin equation (cf. Eq. (pS])) 

U = fa{u) + ^I^)m. (50) 

where force faiu) and noise intensity laiu) are as follows: 

fa{u) = -M + SoU^daiu), 

h{u) =Iu+ (/. + Isu"") dl{u), da{u) = {1 + u^'^y' . (51) 
Corresponding distribution (cf. Eqs. (0), (pO|)) 

Pa{u) = ^e^v{~Ea{u)} (52) 

with partition function Z is determined by effective potential 
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Extremum points of this distribution are determined by equation 



2a/su2" + (l + 



S^-u^-''{i+u^'')] =2a{Is-2Q, 



which gives the boundary of the flux state 

Is = 2/,;, 

related to m = 0. Critical values of state parameters are fixed by condition 
arrive at additional equation 



ld5o I 



(54) 



(55) 



oo to 



u 



2{ 



u 



2a 



u 



-2a 



2als. 



(56) 



- (|56D generalize the simple equalities (|3T|) , (|42D and (146|) related to the 



Expressions 
case a = 1. 

Above expressions show that qualitative results of Section |T| obtained for the particular 
case a = 1 are kept invariant at passage to general case < a < 1. Indeed, the most essential 
difference is observed for noiseless case, namely the steady-state velocity u becomes nonzero 
within whole interval of the driven slope 5*0 (see Fig. 5). Increase of the vertical velocity noise 
J„ causes monotonic w-growth, whereas J5-increase arrives at an effective barrier formation 
near the point m = 0, so that the dependence u{So) becomes non-monotonic in character 
at magnitudes Is above the straight line (|55D (see Figs. 6). Here, by analogy with noiseless 
case (see Fig. 1), lower branches of curves correspond to unstable magnitudes of the order 
parameter, the uppers - to stable ones. According to Fig. 7, the domain, where avalanches 
can not be created, is located near an intermediate magnitudes of the state parameters S'o, 
J„, Is- The phase diagram related to the avalanche formation reveals the same form as for 
the simplest case a = 1, but the straight line (|42|) shifts abruptly to (|55D with escaping 
the point a = 1 (compare Fig. 8 with Fig. 4). According to Fig. 9, increase of the vertical 
velocity noise ly increases the domain of the avalanche formation. 



V. SIZE DISTRIBUTION IN SELF-SIMILAR ENSEMBLE OF AVALANCHES 

Contrary to the previous, when the process of a single avalanche formation has been con- 
sidered, further we aim to study analytically self-similar size distribution (|^) over avalanche 



ensemble. This means that, along line of Section |T|, we will account for noises of complete 



set of freedom degrees, on the one hand, and the fractional feedback type of introduced in 



Section |^ on the other one. Thereby, the Lorenz system unified in above manner is the 
basis of our examination. However, instead of visible geometric-and-mechanic characteristics 
of 'real' sand pile, the system under consideration is parametrized now by a set of pseudo- 
thermodynamical values, which describes the avalanche ensemble in a spirit of the famous 



Edwards paradigm [^, generalized to nonstationary system. In this line, we study 
time dependencies of the avalanche size, nonextensive complexity and nonconserved energy 
of the moving grains. Within the framework of usual synergetic approach, these degrees of 
freedom play roles of order parameter, conjugate field and control parameter, respectively. 
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It is very important that using slaving principle of synergetics and fractional nature of 
the system feedback are shown above to stipulate multiplicative character of noise. It will be 
appeared below, this causes a nonextensivity of applied thermodynamical scheme, so that 
we have to use g-weighted averages instead of usual ones. So, energy of moving sand grains 
is defined by expression 

C^EOp?, (57) 

i 

where Pi is a probability to move grain i with energy ^j, g 7^ 1 is a positive parameter 
being measure of the system nonextensivity determined below. Respectively, nonextensive 
complexity of moving sand grains is an analog of Tsallis entropy ||2^ to be determined as 
follows: 

= -^f^- (58) 



Three-parameter set of standard synergetic scheme is completed by avalanche size s. 

Following above elaborated line, we postulate that self-consistent behavior of the system 
under consideration is presented adequately by set of above quantities governed by the 
Lorenz-type equations (cf. Eqs. (^)) 



TsS = -S + OsEg + ^(t), 

rsS, = -S, + asS^/'C, + (59) 

rcC, = (C° - Q - a^s^/'j:, + JT^ m- 



Here r^, rs, note relaxation times of corresponding values, a^, ae, related feedback 

parameters, Ig, Is, /(■ are respective noise intensities, r is a positive exponent and C° is exter- 
nally driven energy of sand motion. The distinguishing feature of the first of these equations 
is that in noiseless case genuine characteristics s, are connected in linear manner. On 
the other hand, connection of thermodynamic- type values (q, Eg with the avalanche size s 
is stated by two last equations (|59|) to be nonlinear in nature. Physically, this means linear 
relation between complexity and avalanche size near steady-states. Running away these 
derives to negative feedback of the avalanche size with the complexity that arrives at the 
energy decrease - in accordance with Le Chatelier principle, and positive feedback between 
the avalanche size and the energy that brings the complexity increase to be reason of the 
avalanche ensemble self-organization. 

To analyze the system (p9D, it is convenient to measure the time t in unit r^, as well as 
the values s, Eg, (g and Is, /e, are related to the scales: 



1 ^„„ 1 / N 1 



/f = (aEflc)"^, = a;^(asac)~-, /f = aj'^a^^^^^^ a^^^ ^\ (60) 
Then, the modified Lorenz system (^) takes the simple form 

s = -s + Sg + /^e(i), 

^t, = ~^g + s^/\g + ^m, (61) 
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where relaxation time ratios are introduced: 



^ = r^/r,, 9 = tJts. (62) 



It is worthwhile to notice that system (plf ) is passed to the form of Eqs. (^) if the values s, 
Eg, (^g, r/2, 6 are replaced by m, f , S', a, e, 5, respectively. 

It is well-known that complete set of SOC systems can be reduced to one of two fam- 
ilies [|r^: systems with deterministic dynamics extremely driven by random environment 
(growing interface models, Bak-Sneppen evolution model and so one) and the stochastic 
dynamics family (models of earthquakes, forest-fire et cetera). [] Remarkable peculiarity of 
obtained system (|6T|) is a possibility to present in natural manner both mentioned families. 
So, the former is related to noiseless case, when 7^, J^, = but magnitude of the energy 
relaxation time is above than the ones for complexity and avalanche size (r^ > r^, r^); on the 
other hand, a parameter of environment drive C° has to take magnitudes above the critical 
one Cc = 1 [0- In such a case, the system (^) describes strange attractor that may 



represent behavior of the first type SOC systems. Proper stochastic behavior is relevant 
to nonzeroth noise intensities Is,I^,I( 7^ that make possible the SOC regime appearance 
even in absence of driven affect {(^ = 0). 

Taking into account that problem of Lorenz strange attractor is well known P? 



we 



will restrict ourselves further with treatment of the stochastic system, where the adiabatic 
conditions i),6 <^ 1 are applicable. Then, two last equations of the system ( pT] ) arrive at 
dependencies type of Eqs. (p3D 



Sg(t) = Sg + Sge(t), Ut) = C, + C,at), (63) 
where deterministic and fiuctuational components are determined as follows (cf. Eqs. (|2^)) 



(64) 

4 = C° dr{s), Cg = ^Ic + I^s-d^is), d^s) = {l + s^-'. 

So, due to the slaving principle of synergetics, initially adiabatic noises of the complexity 
and energy are transformed to multiplicative form. On the other hand, relation between the 
complexity and energy 



= VC.(C°-C,), (65) 
following from the dependencies ([B^), arrives at expression 

(66) 

for effective temperature T = dQq/dTjq. As is depicted in Fig. 10a, so defined temperature 
increases monotonically with energy growth from magnitude T = at = to infinity 




^Generally, we face here with much more complicated problem, see Ref. |31|. 
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at the point (q = C°/2- Here, the temperature T breaks abruptly to negative infinity and 
then increases monotonically again to initial magnitude T = at (q = This means that 
inside domain < Q < the avalanche system is dissipative to behave in usual manner; 
contrary, within domain (^/2 < Cq < C° self-organization process evolves, so that an energy 
increase derives to complexity decrease, in accordance with negative value of temperature. 
At steady state, where avalanche has stationary size sq = y/C'^ — 1, the temperature takes 
the stationary magnitude 



Jr^ - 1 

being negative within supercritical domain 1 < C° < 2. According to Fig. 10b, the stationary 
temperature Tq decreases monotonically with the driven energy from the zeroth magnitude 
at C'' = 1 to negative infinity at = 2. 

Presented self-organization regime addresses to externally driven systems, which are rel- 
evant to the usual phase transition but not to the SOC itself. To study the latter within 
above consideration, let us combine Eqs. (|6^), (Q) with the first of equations (0) along the 
line, which has been used above for obtaining Langevin equation (p5l). So, by analogy with 



Section we arrive at stochastic equation (|50|), where effective force and noise intensity 
are given by Eqs. (|5l|) with accuracy to the replacements mentioned after Eqs. (|6^) : the 
quantities s, S,, (q, t/2 have to be taken instead of u, v, S, a. Then, all results obtained in 
Section |I^ can be used immediately. Particularly, it is appeared that infiuence of random 
scattering of avalanche size is non-essential, whereas energy and complexity noises affect cru- 
cially. Related picture is reflected by Fig. 8 taken in plane — Is formed by corresponding 
noise intensities of avalanche ensemble. By this, mixed domain A+N respected to intermit- 
tency regime is bounded by the straight line (|55D and bell-shaped curve type of Eqs. (^6]). 
According to Fig. 9, where exponent a has to be replaced by t/2, random scattering growth 
of the complexity extends the SOC domain along axis of exponents r. 

We are in position now to consider the avalanche size distribution on the basis of equa- 
tions ( |5TD - (1531). At arbitrary relations of noise intensities, general expressions are as 
follows: 



n^) = f(Sexp|/|g}d.'}; 



/(s) = -s + CV/2rf,(s), (68) 



I{s) = h + (Is + I^s-) dUs), d^s) = {1 + s^r' 



At SOC regime, the driven energy C° = and the distribution ( |68D behaves as depicted 
in Fig. 11 for different noise intensities of both energy and complexity. It is seen that the 
power-law dependence inherent in SOC regime is observed only within limits s ^ 1 and 
<^ /( (at condition = 0). Here, the distribution (|68|) is reduced to the canonical 
form (0), where second multiplier takes the form 

P(.) = ^exp|-/,-Y^^d.j, dAs)^{l + sT'. (69) 
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It is easily to see that deviation of this multipher off a constant value is estimated with term 
~ s^""^, whose contribution increases remarkably with r-decrease and growth of avalanche 
size to extremely large magnitudes s ~ 1, i. e., with escaping SOC domain. This confirmed 
by Fig. 12, where deviation St of the slope of dependence P{s), Eqs. (|68|) within linear 
domain from the theory parameter r is depicted as function of the parameter r itself. It 
is seen, in accordance with above estimation, that the deviation 6t takes maximal value 
6t < 10~^r at non-essential magnitudes r < 1 or with noise intensity growth to enormous 
values Is ~ 10'^. 



VI. DISCUSSION 



Remarkable peculiarity of expression (|69|) is that, within limits s -C 1, /s, /s <^ -^c 
inherent in SOC regime, it can be rewritten in the form 



dzH 



exp 



r(2 



(70) 



expressed in terms of standard Gamma-function T{x) and fractional integral 1"^/ of order 
2 — r, whose definition ( |A.1| ) is given in Appendix (for details, see Refs. [0, [0). On the 



other hand, it is well-known |^ that such type expressions appear as a solution of fractional 
Fokker-Planck equation 



V^V{s,t)=V^AsV{s,t) + 



Tiw) 



(71) 



where fractional derivative defined by Eq. ( [A.2| ) is used to be inverted to the fractional 
integral (|A.1|) . Multiplying equation ([7T|) by term s^^, for a a-average 



|s| = (s")^, (s") = y s"P(s,t)ds, a>0 



(72) 



one obtains at a = 2vu 



t, z 



2zu 



(73) 



where z is dynamical exponent. This relation corresponds to large time limit, where only 
diffusional contribution is essential. Combining expressions ([70|), ([7^) and ([A. 1|) arrives at 
relations 2 — t = w = zuj/2^ which yields 



ZUJ 



Comparing this equation with the second of known relations (|1]), one obtains 



IxJZ 



D' 



(74) 



(75) 



Then, mean-field magnitudes a; = 1 and D = A are related to dynamical exponent z = 1 
that, in accordance with definition (^), addresses to unusual ballistic limit of SOC regime. 
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On the other hand, the fractional Fokker-Planck equation ( [fTD arrives at usual diffusional 
regime with z = 2 only in artificial case, when time-derivative exponent is supposed to be 
u = 1/2. 

Obvious reason for such discrepancy is non-consistent application of usual field relations 
(|I|) to Lorenz system (0). In this system, the role of different space directions is played 
by stochastic degrees of freedom s, Eg and (g, whose number is = 3. However, stochastic 
process evolves for every of these variables in plane spanned by given variable itself and 
conjugated momentum. Moreover, multiplicative character of noise, which is determined by 
exponent a in expressions (^Tj), reduces fractal dimension of every plane to the value 2(1 — a) 
p8| . Thus, resulting fractal dimension of phase space, where stochastic system evolves, is 
as follows: 



D = 2n(l - a) 



(76) 



where one has n = 3 for used Lorenz system. Insertion of this dimension into expression 
(fTSP arrives at the connection uz = 2, which in contrast to above obtained relation uz = 1 
is correct in the simplest case u = 1, z = 2 [the latter is relevant to single stochastic freedom 
degree (n = 1) with additive noise (a = 0)]. In general case, equations ( |7^ - ([7^) yield final 
result 



r 



1 - 



1 



2n{l - a) 



(77) 



Respective dependencies are depicted in Figs.l3a,b to show that exponent r increases mono- 
tonically from minimum magnitude r = 1 at critical number (1 — a)~^ to upper value r = 2 
in limit n ^ oo; thereby, a-growth shifts dependence r(n) to large magnitudes n, i. e., 
decreases exponent r. 

It is easily to see that relation ( |77| ) reproduces known results of different approaches 



for dimension D (see Ref. [0]). In the case related to mean- field theory, one has r = 3/2 
and equation (|77D expresses number of self-consistent stochastic equations being in need of 
treating SOC behavior as function of exponent of the corresponding multiplicative noise: 



n 



1-a 



(78) 



Thus, in accordance with Fig. 13c, self-consistent mean-field treatment becomes possible 
if the number of relevant equations is more than minimum magnitude = 2. Approaches 



0' El 



1^ represent examples of such considerations, where noise is supposed to be 
additive in character (a = 0). Switching multiplicative noise arrives at a-growth and non- 
contradicting representation of SOC demands increasing number of self-consistent equations: 
for example, within the field scheme |20| related to directed percolation (a = 1/2), mean-field 
approximation is applicable for dimensions more the critical magnitude dc = 4; used here 

3) addresses to multiplicative noise characterized 



and in Refs. IT^, ||T8[ Lorenz scheme (n 
by the exponent a = 1/3 (see below). 

Let us focus now on relation of above exponents to nonextensivity parameter q addressed 
to Tsallis definitions 



24]. Relevant kinetic behavior could be described by non- 



linear Fokker-Planck equation 
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(79) 



where measure units are chosen so that effective diffusion coefficient disappears, 00 > 0, 
q > are relevant exponents [Q, IQ. Supposing here normahzed distribution function in 
self-similar form type of Eq. (|^) 



P{s,t) = s^'^V{x); Sc = Sc{t), x = s/sc 



(80) 



we obtain 



si ~ r, V"-' ~ X. (81) 

On the other hand, we could use the fractional Fokker-Planck equation type of Eq. (ffTI): 

V'^P{s,t)=V'_^^P{s,t). (82) 
Inserting here the solution (BO), one finds dependencies 



2ro jo; 'T-)2ro— 1 



whose comparison with Eqs. (|5I|) yields 



(83) 



(84) 



Because of the average |s| in Eq. i\7\ 
systems, relevant dependencies (|73| 



is reduced to the scale Sc in the case of self-similar 
ID and m) get 



q = zuj. 



(85) 



Combining this equality with Eqs. (|75D, ( [7B| ) arrives at the resulting expression for nonex- 
tensivity parameter of the system under consideration: 



n{l — a) 



(86) 



Maximum magnitude q = 2/n is related to systems with additive noise (a = 0), which 
is relevant to mean-field picture at number of governing equations < 2 (n = 1, really). 
Switching multiplicative noise with increasing exponent a > arrives at g-growth and self- 
organizing system becomes nonextensive in character {q > 1) at the boundary ( [78| ) of the 
mean-field applicability domain. In accordance with above estimation, fractional Lorenz 
system is nonextensive essentially if the exponent a > 1/3. 

It is worthwhile to remember that we have addressed above the superdiffusion process 
only to be related to Levy flights at discrete time instants with arbitrary displacements, 
including infinite ones Being related to the Fokker-Planck equation (0), such type 

processes are characterized by exponents 00 = 1 and w < 1, the first of which is constant, 
whereas the second one characterizes fractal time-sequance of the Levy flights to arrive at 
the dynamical exponent z = 2zu/u < 2 (see the last of equations ([73|) ). To be dependent 
on microscopic conditions, the probability distribution of displacements x reads 
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(87) 



that is characterized by the fractal dimension D and microscopic step exponent 7. It is 
appeared that in the case of rare events, when 7 < 2, the dynamical exponent z is reduced 



to the microscopic step exponent (2; = 7 < 2), whereas at 7 > 2 one has 2: = 2 |36 . 

In the opposite case of subdiffusion process, a microscopical ingredient is random Levy 
walks instead of the discrete Levy flights. This process is involved continuously in course of 
the time over discrete placed traps, so that the exponents u < 1 addresses fractal properties 
of this space to be dependent on microscopic conditions. This properties produce trans- 
formation of the usual Boltzmann-Gibbs statistics in nonextensive manner [|^. Along this 



line, the subdiffusion process is presented by the Tsallis-type distribution [3/ 



p(x) oc [1 + I3{q - \)x^\ 13 = const > 0, (88) 

where deviation g — 1 of the nonextensive parameter is caused by fractal nature of the system 
phase space to be connected with the step exponent 7 as follows: 

Formal advantage of the distribution (^8D is that corresponding g-weighted average 

(x^), ^ / xV(x)d^x, (90) 



where the integrand varies as x (^+'^), is converges for arbitrary step exponents 7 > 0. As a 
result , the law of motion of the random Levy walker is given by 

^^^^^^ ' '^ = il-(g-l)f at 7>2, g > 1. 



In contrast to Eq. (p3|), where the exponents w < 1, uo = 1 are relevant to the super diffusion, 
here one has inverted relations w = 1, uo < 1. Thus, in accordance with subdiffusion nature, 
the last equahty ( [73| ) yields the dynamical exponent z > 2. 



In general case -07,0; 7^ 1, inserting Eqs. (|9TD into the relation (|85|) arrives at the result 



if,„ at 1 < g < gD, 



i-f (5-1) . - ^g2) 

^ at gz3<g<2, ^ ' 

where boundary value of the nonextensivity parameter is introduced: 

A + D 

to - J^. (93) 

To avoid a mistake, let us focus that in contrast to equalities ([1|), (|7^, (|75D, which could be 
addressed to both the real phase space and the configurational one (the latter is spanned by 
variables of governing equations), above obtained relations ( PT| ) - ( PBD are relevant to the 
real phase space only. This is reflected by addressing the fractal dimension D to the only 
real coordinate space in the former case, whereas in the latter it is reduced to the effective 
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value (ffGl). Along the principle line of our treatment, central role is played by the relation 
(76) since, by analogy with studying renormalization group, we have considered mainly the 
properties of the configurational space but not real diffusion process. 

Finally, let us address the question: Why the Lorenz system is used but not Rossler 
one or another? The reason can be recognized within supersymmetry field approach, whose 
using shows that the Lorenz system could be generated by the Langevin equation for order 
parameter to be relevant to standard stochastic system . On the other hand, one can see 
within microscopic consideration that Lorenz system addresses to the simplest Hamiltonian 
of a bozon-fermion system [^. To convince in this statement, let us consider such system 
with interaction w. Bozons are described by creation and annihilation operators hf , hi, 
satisfying the usual commutation relation: = ^im, where /,m are the site numbers. 

The two-level Fermion subsystem is represented by operators af^^aia-, a = 1, 2, for which 
the anti-commutation relation {a;^, a^^} = SimSafs is fulfilled. The occupation numbers 6j^6k 
determine the Bozon distribution within k-representation that corresponds to the Fourier 
transform over lattice sites /. To represent the Fermi subsystem we should introduce the 
operator di = af^ai2 determining the polarization with respect to the saturation over levels 
a = 1, 2, as well as the occupation numbers nia = a^a/^. As a result, behavior of the 
system is defined by Dicke Hamiltonian 

^ = E { (^i^ki + ^2^k2) + oj^b+b^ + '-w{b+d^ - 46k) } , (94) 

where Ei^2 are energies of the Fermi levels, cuk is the Bozon dispersion law, the imaginary 
unit before the interaction w reflects the Hermitian property and the Planck constant is 

n = i. 

The Heisenberg equations of motion corresponding to Hamiltonian (p4D have the form 



w 

6k = -ic^k^k + —dk] (95) 

w 

c?k = -iAdk + — 6k(?^k2 - ^ki), A = E2-Ei; (96) 
w w 

^^kl = ^l^kC^k + 4^k), ?^k2 = - — (6kC?k + rfk^k)- (97) 

At condition of resonance, the first terms in the right hand sides of equations (0), ( P7D con- 
taining frequencies Wk and A may be suppressed by introducing the multipliers exp(— iwk^) 
and exp(— iAt) for the time dependencies 6k (^), d\f^{t), respectively. On the other hand, to 
take into account the dissipation, these frequencies acquire imaginary additions — i/r^;, — I/t^ 
characterized by relaxation times r^, Ty (here the conditions Im c^k < 0, Im A < reflect 
the causality principle). As a result, equations (pUf), ( |97D get the dissipative terms —b^/T^, 
—dk/Ty, where is the relaxation time of Bozon distribution and Ty is the Fermion po- 
larization time. Obviously, one can suppose that the dissipation influences onto the Fermi 
levels occupancies nka(i) also. However, since the stationary values n^^, 7^ (in case of 
external drive > ^ki) ^^le dissipative terms in Eqs. (0) take much complicated form: 
— {niia ~ ^\^a)/'TSi where ts is the relaxation time of the Fermion distribution over levels 
a = 1, 2. 

Now, let us introduce the macroscopic quantities 
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= (bt) = (&k), V]^ = (rfk) = (c?k), 
5k = K2-nki), 5^=K2-<i), (98) 

where the angular brackets mean thermodynamical averaging. Then, neglecting the correla- 
tion in distribution of particles over quantum states and omitting dependence on the wave 
vector k, the Heisenberg equations ( p^ ) - (pTf), being complemented by dissipative terms, 
result in the initial Lorenz system (^, (^, (^), whose parameters take magnitudes a = 1, 
Tx = 2/w, D = {2w)~^. Respectively, the dimensionless variables in the system (|TD|) - (|T2p 
are as follows: u = w^TyTsY^'^x, v = w{TyTsy^'^y and S = {wTy/2)y'. 

As a result, we arrive at the conclusion that Lorenz system is relevant microscopically 



to the simplest bozon-fermion system defined by Dicke Hamiltonian (|9^). At first glance, it 
could have been shown that one can write corresponding expression for macroscopic (effec- 
tive) Hamiltonian being a synergetic potential, whose dependence on the freedom degrees u, 
V, S could have generated the Lorenz system. But such dependence is appeared to be absent 
because of the effective Hamiltonian has to take into account quite different commutation 
rules related to different freedom degrees. Obvious advantage of above mentioned super- 
symmetry theory and stated microscopic approach is that we have explicit possibility to 
account such difference. Generally, this situation is relevant to known problem of description 
of systems with intermediate statistics (see Ref. 
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VIII. APPENDIX 

Here the basic properties of fractional integral and derivative, as well as Jackson deriva- 
tive are quoted for convenience. The integral of fractional order w is defined by equality 



where f{x) is arbitrary function, T{x) is standard Gamma-function. To be inverted to the 
fractional integral, relevant derivative = of order > is determined as follows: 



At < -07 < 1 it is convenient to use equation 

X 

zu r 
T{l-zu)J {x-x') 



™-^7^/#^d.'. (A^3) 
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where we take into account known equality xr{x) = r{x + 1) for x = —w. 

Let us introduce now Jackson g-derivative, whose advantage for analyzing self-similar 
system is that this derivative determines the rate of a function j[x) variation with respect 
to dilatation g 7^ 1, but not to the shift dx — > 0, as in usual case g = 1. According to such 
definition, Jackson g-derivative reads: 



g > 0. 



For important case of homogeneous function, being subjected to condition 

f{qx) = q'^fix), 



(A.4) 



(A.5) 



where g > is a dilatation parameter and a > is an exponent, the Jackson g-derivative is 
reduced to Jackson g-number: 



'^qfix) = [a]J{x), [a\ 



a'' 



(A.6) 



It is easily to see that the value [a]q — > a in the limit g ^ 1 and varies as g" ^ at g ^ cxd. 
On the other hand, the Tsallis g-logarithmic function 



lHqX 



g-1 



(A.7) 



can be represented in the form of the Jackson g-number ( [A.6|) with the index a - 
l)(lna:/ Ing). Accompanied by Eqs. ( |A.6|) this relation and obvious equality 

\ng{xy) = \nqX + lnqy + (g - 1) {\nq x) (\ng y) 

lead to important rule for the Jackson derivative: 

Vq [f{x)gix)] = [Vqfix)] g{x) + f{x) \Vqg{x)] + (g - 1) \VJ{x)] \Vqg{x)] . 

Generalizing Eqs. ( [A.3|) , ( [A.4| ), let us introduce finally a fractional Tu-derivative: 



(A.8) 



(A.9) 



wx 



Til-vj) 



f{x) - fjqx) 
(1 - g)i+^ 



dg, <w < 1. 



(A.IO) 



In the case of self-similar system, the function /(x) is homogeneous to be subjected to 
condition ([A.5|) . Then, the definition (|A.1CI|) is simplified: 



V^f{x) = {a}^x-^f{x), 0<w<l 
to be reduced to fractional cc-number 



(A.ll) 



r 1-w) 



q^ 



l+VJ 



dg, < a, < < 1. 



(A.12) 
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Being the combination of F-functions: 

r(a + i) 



1 



T(a + l-uj) T(l 



VJ] 



(A.13) 



this number increases monotonically with growth of both exponents a and w taking zeroth 
magnitude on the axes -07 = 0, a = and characteristic values \ \\ ^ ~ ^ ~ — 0-322, 

= 1. Such behavior is characterized by the particular dependencies 



a 

Til 



a] 



1 



r(l-a) 



at w 
at w 



a, 



(A. 14) 



r(2-ro) 



at a = 1. 



Thus, if g-number ( |A.6|) related to Jackson derivative (|A.4|) tends to exponent a in the limit 
g — > 1, -07- number ( |A.12| ) corresponding to fractional integral (|A.10|) is reduced to factor a 
at Ct7 = 1. 
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FIGURE CAPTIONS 



Fig. 1. The S'o- dependencies of (a) the velocities Ue, u"^, and (b) the equihbrium slope 
Se- The arrows indicate the hysteresis loop. 

Fig. 2. Phase portraits in the v — u plane at m = 1, mq = 0.1, S'o = 1.25S'c: (a) e = 10^^; 
(b) e = 1; (c) e = 10^. 

Fig. 3. Phase diagrams at fixed values ly-. (a) 1^ = 0; (b) = 1; (c) = 2. Curves 1 
and 2 define the boundary of stability of avalanche (A) and non-avalanche (N) phases. 

Fig. 4. Phase diagram for system with Sq = and Ig, 1^ ^ {D - disordering point; T 
- tricritical point; C - critical point). 

Fig. 5. The So-dependence of the steady-state velocity u at a = 0, 0.5, 0.7, 0.9, 1.0 
from top to bottom. 

Fig. 6. The S'o-dependence of the steady-state velocity u: (a) at a = 0.75, = 1 
(curves 1-4 address to Is = 1, 2, 3, 5); (b) at = 1, Is = 5 (curves 1-4 address to 
a = 0.25, 0.5, 0.75, 1.0. 

Fig. 7. Three-dimensional phase diagram (the non-avalanche domain is located under 
surface) . 

Fig. 8. Phase diagram for system with Sq = and Is, Iv ^ at a = 0.5, 0.75, 1.0 
(dotted, solid and dashed curves, respectively). Diamonds address to curves 1 - 4 in Fig. 11. 

Fig. 9. Phase diagram in the Is — a plane at J„ = 2, 3, 4, 5, 6 from bottom to top (the 
non-avalanche domain is located inside the curves). 

Fig. 10. The energy dependences of the avalanche ensemble temperatures: (a) nonsta- 
tionary magnitude T versus ratio Cg/C^i (b) stationary temperature Tq versus 

Fig. 11. Distribution function ( |68D at r = 1.5 and regimes pointed out by diamonds in 
Fig. 8: 1) h = 0, Is = 50 (SOC); 2) I, = 0.5, Is = 30 (A+N); 3) = 1, /cj = 5 (N); 4) 
4 = 2, Is = 0.5 (A). 

Fig. 12. Deviation 6t of the linear slope of curve 1 depicted in Fig. 11 from parameter 
r versus the exponent r itself {Is = 10, 50, 10^ from bottom to top). 

Fig. 13. Dependences of exponent r: (a) on equations number n {a = 0, |, |, | from top 
to bottom); (b) on exponent a {n = 2,3,4,6,8, 10 from bottom to top); (c) phase diagram 
for mean-field and nonextensivity domains. 
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